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The observation of a narrow structure near the J /ip </> threshold in exclusive iJ+ — > 
J/tlxpK'^ decays produced in pp collisions at ^ = 1.96 TeV is reported. A signal of 
19 it 6(stat) it 3(syst) events, with statistical significance of 5.0 standard deviations, is 
seen in a data sample corresponding to an integrated luminosity of 6.0 fb~^, collected by 
the CDF II detector. The mass and natural width of the structure are determined to be 
m = 4143.41^'^ (stat) it 0.6(syst) MeV/c^ and F = IS.SlgVC^t*^*) ± 2.5(syst) MeV/c^, 
consistent with the previous measurements reported as evidence of the y(4140). 

PACS numbers: 14.40.Rt, 13.25.Gv 

1. Introduction 

In recent years, several states with charmonium-like decays were discovered, which 
do not fit properly into the established charmonium picture. These states, called 
X, Y, Z, are candidates for exotic mesons beyond the conventional quark-antiquark 
model (qq). Possible interpretations are quark-gluon hybrids {qqg), four-quark states 
{qqqq), molecular states composed of two usual mesons, glueballs, et cetera. 

CDF reported an evidence for a narrow near-threshold structure in the J/ipct) 
mass spectrum, called y(4140), using exclusive J/ip(j)K^ decays.^As there 

was no signal seen by Belle in a subsequent searc h,^ it is important to investigate 
with a larger CDF data sample. Consisting of two vector mesons (positive C-parity), 
the final state J/ip(j) is a good channel to search for an exotic meson. The observed 
excess would be the first charmonium-like structure decaying into two heavy quarko- 
nium states (cc and ss). Because the mass in this channel is high enough for open 
charm decays, an explanation as charmonium state is very unlikely for a narrow 
structure. A search near the J/ip4> threshold is also motivated by the closeness of 
the y(3930) to the J/ipuj threshold. Some of the possible exotic explanations for 
the F(4140) are discussed in Ref. ISl 
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2. Candidate Selection 

We report on an update of the search for structures in the J/tp4' system produced in 
exclusive — J /ijjclyK^ decays.^ The employed dataset was collected by the CDF 
II detector at the Tevatron and corresponds to an integrated luminosity of 6.0 fb~^. 
It was accumulated using a dedicated dimuon trigger which requires a /i'^/x" pair 
with a mass of 2.7 < m(/i+^~) < 4.0 GeV/c^. Due to trigger prescales at increasing 
instantaneous luminosities one cannot expect a linear increase of the sample size 
compared to the previous analysiJJ with 2.7 fb^^. 

In order to build 5+ — )■ J/ipcjiK'^ candidates, first J/ip ^ l^'^f^^ and (j> — > 
candidates are reconstructed which are then combined with an additional 
charged track with kaon mass hypothesis. Thereby, the reconstructed J/ip and (jj 
masses are required to lie within 50MeV/c^ (J/i^) respective 7MeV/c^ {(f>) of the 
corresponding world average values. The combinatorial background can be reduced 
significantly with a higher threshold of the decay length in the transverse plane 
{Lxy{B^)) to exploit the long B meson lifetime. In addition, a kaon identification 
quantity can be used to obtain a further background reduction by separating the 
final state kaons from the dominant pion background. For that purpose, the in- 
formation about the ionization energy loss dE/dx in the drift chamber and the 
information from the Time-of-Flight detector are combined in a log- likelihood ratio 
LLRKaon- The cuts Lxy{B~^) > 500 /Ltm and LLi?Kaon > 0.2 are chosen by opti- 
mizing the quantity S/y/S + B, where 5* and B are the numbers of B~^ signal and 
background events, respectively. After both requirements, a background reduction 
factor of approximately four orders of magnitude is accomplished. 

Figure [Ija) shows the resulting J/^xf'K'^ mass spectrum. A fit with a Gaussian 
signal and a linear background function yields 115 ± 12 signal events, corresponding 
to a 53% increase over the previous analysis. For the examination of the J/ip(j) spec- 
trum, only candidates within ±3(j (±17.7MeV/c^) around the nominal B^ mass are 
selected. Furthermore, sideband events within [—9, -~6]a and [-1-6, +9]a around the 
nominal _B+ mass are used to model the combinatorial background in the J/ip(j) spec- 
trum. Figure [ijb) shows the mass difference AM = m{fi^ fi~ K'^ K") — m{fi'^fi~) 
distributions of the resulting J/4'4' candidates. Whereas the y(4140) can be seen as 
narrow near-threshold excess in the B mass window, no such evidence is found from 
the B mass sidebands. Figure [2l[a) shows the Dalitz plot of the candidates from the 
B mass window and figure [JJb) the B sideband-subtracted K^K~ mass spectrum 
without mass window requirement. The fit function is a P-wave relativistic Breit- 
Wigner convolved with a Gaussian to account for the detector resolution. As there 
is no significant background contribution, the B^ — > J/ipK~^K~K^ final state is 
well described as J/ip(j)K'^. A comparison between the J/tpcj) mass difference dis- 
tributions of the dataset used for the updated analysis described in this write-up, 
corresponding to an integrated luminosity of 6.0 fb~ , and the one employed for the 
published y(4140) measurement"'^, corresponding to 2.7 fb""'^, can be found in figure 
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Fig. 1. (a) J/ip(f)K^ mass distribution with a fit to the data represented by the soUd blue Une. 
The vertical dashed black and red lines indicate the mass window and sidebands described 
in the text, (b) Mass difference distributions of the resulting J/tp<j> candidates from the mass 
window (black histogram) and sidebands (red histogram). 




Fig. 2. (a) Dalitz plot of the final state J/i/'(/)ii'+ in the B+ mass window. The boundary shows 
the kinematically allowed region, (b) K~ mass distribution with the fitted function (solid blue 
line) as described in the text. 
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Fig. 3. Comparison between the J/tl)<t> mass difference histograms of the dataset used for the 
updated analysis described in this write-up (6.0 fb~^) in black and the one employed for the 
published y(4140) measurement (2.7fb~'^) in dashed red. 

3. Fits to Data and Significance Determination 

Figure HJa) shows the AAf distribution from the B mass window excluding events 
with AM > 1.56 GeV/c^ in order to avoid combinatorial backgrounds from misiden- 
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Fig. 4. (a) J/ip'P mass difference distribution with a fit to the data represented by tfie solid 
red line, (b) —2ln{Co/ Cmax) distribution for 84 million simulation trials. The vertical red line 
indicates the value obtained in data. 



tified — > ip{2S)(l) — ?> {J/ip tt^ tt^) (j> decays. An unbinned maximum likelihood 
fit is performed, where the enhancement is described by the convolution of an S*- 
wave relativistic Breit-Wigner function with a Gaussian resolution of 1.7MeV/c^ 
obtained from Monte Carlo simulations, and the background is modeled by three- 
body phase space. Including systematic uncertainties, which are estimated by vary- 
ing the fit model, the fit yields 19±6(stat)±3(syst) signal events. It returns a mass of 
m = 4143.4it3 o(stat)±0.6(syst) MeV/c^ after including the world average J/ijj mass 
and a decay width of F = 15.3tg°;^^(stat) ± 2.5(syst) MeV/c^, both consistent with 
the values from the published measuremenlfH. The observed width, which is much 
larger than the resolution, suggests a strong decay for the y(4140). Furthermore, 
the relative branching fraction to the nonresonant — > Jf-tjj (p decay is mea- 

sured as ''^''^"^Xb?-^// = 0.149 ± 0.039(stat) ± 0.024(syst), where 
the relative efficiency is determined to be 1.1, using an S'-wave relativistic Breit- 
Wigner function with mean and width values determined from data to represent 
the y(4140) structure and three-body phase space kinematics for the nonresonant 
B+ ^- J/i}(l)K+ decay 

In order to estimate the probability of a creation of such a signal due to back- 
ground fluctuations, a large number of three-body phase space decays are per- 
formed and the number of trials which produce a signal with a log-likelihood ratio 
— 2h\ [Cq / Cmax) of the null hypothesis fit and the signal hypothesis fit larger than 
the value measured in data are counted (see figure [UJb)). Thereby, the mass can be 
anywhere in the considered AM window and the width has to be larger than the 
detector resolution and smaller than 120 MeV/c^. This procedure leads to a p- value 
of 2.3 • 10"^, corresponding to a significance of the enhancement of 5. Oct. 

In figure [UJa), an additional excess above the background appears at a mass of 
approximately 1.18GeV/c^. With the parameters of the F(4140) fixed to the val- 
ues obtained from the fit described above, another unbinned maximum likelihood 
fit assuming two structures and the same background model as before is performed. 
Thereby, the additional enhancement is described by the convolution of an S'-wave 
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Fig. 5. (a) J l^4> mass difference distribution with a fit to the data, where an additionai signai 
structure is inciuded which is located about one pion mass higher than the y(4140). (b) Distribu- 
tion of the log-likelihood ratio of the fit assuming only the y(4140) and the fit assuming two signal 
structures for the simulation trials. The vertical red line indicates the value obtained in data. 

relativistic Breit-Wigner function with a Gaussian resolution of 3.0MeV/c^ ob- 
tained from Monte Carlo simulations. The measured data distribution together with 
the fit projection can be found in[5l[a). The fit returns a yield of 22 ±8 signal events, 
a mass, after including the world average J/^ mass, of m = 4274.4j;*;7MeV/c2 and 
a decay width of T = 32.3ti5 3 MeV/c^. Just like in the y(4140) case, the statis- 
tical significance of the additional excess is determined by simulations, where the 
log-likelihood ratio of the fit assuming only the F(4140) and the fit assuming two 
signal structures is calculated (see figure[5l[b)). This leads to a p-value of 1.1 • 10"'^, 
corresponding to a significance of 3.1 cr. 

After the first confirmation of Belle's X(3872), including the determination of its 
allowed quantum numbers and the most precise mass measurement,^ CDF keeps 
contributing to the field of exotic X, F, Z states with this recent observation of 
a narrow structure, called F(4140), near the J/i'<t> threshold in exclusive B+ ^■ 
J/%l)(t>K'^ decays. 
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